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ABSTRACT 

We examine the agreement between the observed and theoretical low-mass (< 0.8M©) stellar main sequence 
mass-radius relationship by comparing detached eclipsing binary (DEB) data with a new, large grid of stellar 
evolution models. The new grid allows for a realistic variation in the age and metallicity of the DEB population, 
characteristic of the local galactic neighborhood. Overall, our models do a reasonable job of reproducing the 
observational data. A large majority of the models match the observed stellar radii to within 4%, with a mean 
absolute error of 2.3%. These results represent a factor of two improvement compared to previous examinations 
of the low-mass mass-radius relationship. The improved agreement between models and observations brings 
the radius deviations within the limits imposed by potential starspot-related uncertainties for 92% of the stars 
in our DEB sample. 

Subject headings: binaries: eclipsing — stars: evolution — stars: fundamental parameters — stars: low-mass 
— starspots 



1. INTRODUCTION 

[ The disagreement between the theoretical and observational 
low-mass, main sequence mass-radius (henceforth MR ) rela- 
tionship has been recognized for nearly four decades dHoxid 
119701 119731) . Although, only in the past two decades has 
the disagreement become overwhelmingly apparent with the 
reduction of observational uncertainties (for an excellent re- 
view, see Torres et al. 2010) and the development of sophisti- 
cated low-mass stellar models (Baraffe et al. 1998). 

The primary line of evidence stems from the study of 
detached double-lined eclipsing binaries (hereafter DEBs) 
with additional support garnered by direct measurements 
of stellar radii via in terferometry (e.g., Berger et al. 200^ 
\ I von Braun et al.ll2012 !). These observations routinely suggest 
that stellar evolution models systematically under predict stel- 
lar radii by 5 - 15% and over predict effective temperatures at 
the 3 - 5% level. However, it is presently not clear whether 
the routinely quoted 5-15% disagreement is representative 
of true radius discrepancies or whether there are other factors 
contributing to the derivation of such large radius errors. 
' One such factor derives from the fact that previous studies 
focusing on the comparison between models and observations 
have generally applied a limited sample of isochrones to their 
] data. Largely, these sets are comprised of 1 Gyr and 5 Gyr, 
solar metallicity isochrones. This is predominantly a conse- 
quence of the limited age and metallicity range of currently 
available low-mass stellar models. Age and metallicity ef- 
fects are less important in the low-mass regime, but the strin- 
gent uncertainties quoted by observational efforts preclude 
the use of such a lim ited set of isochrones. For example, 
iBurrows et alJ (1201 ll) discovered non-negligible radius vari- 
ations in brown dwarfs and very-low-mass stars (< O.IM©) 
when allowing for a more comprehensive set of metallicities. 
Isochrones with metallicities spanning a range characteristic 
of the local galactic neighborhood are therefore essential to 
accurately assess the validity of stellar evolution models. 

Furthermore, one must also consider that the population 
of well-characterized DEBs has, until recently, consisted of 

' Neukom Graduate Fellow, Gregoiy.A. Feiden.GR.@Dartmouth.edu 



eight systems. While unlikely, it is not unimaginable that 
those eight systems were more the exception than the rule in 
terms of their lack of consiste ncy with ste l lar evo lution mod- 
els. Since publication of the Torres et al.l (|2010') review, the 
population of well-characterized, low-mass DEBs has more 
than doubled. The availability of this new data allows for a 
more accurate statistical characterization of the agreement (or 
lack of) between the MR relationship defined by models and 
observations. 

When discrepancies are observed, they are typically at- 
tributed to the effects of a large scale magnetic field (e.g., 
Ribasl2006 f: Lopez-Morales 2007; Morales et a l. 2008. 2009ai 
Torres et al. 2010; Kraus et al. 2011) as DEBs are often found 
in tight, short-period orbits with periods under three days. 
Tidal interactions and angular momentum conservation act to 
synchronize the orbital and rotational periods of the compo- 
nents, increasing the rotational velocity of each star in the pro- 
cess. The dynamo mechanism, thought to be responsible for 
generating and sustaining stellar magnetic fields, is amplified 
as a result of the rotational spin-up and enhances the efficiency 
of magnetic field generation within the star. Each component 
in the binary system is then more able to produce and maintain 
a strong, large-scale magnetic field than a comparable single 
field star 

The effects of a large-scale magnetic field are thought to be 
two-fold: convective motions within the star are suppressed 
and the total surface coverage of starspots is increased. In 
both cases, a reduction in the total energy flux across a given 
surface within the star occurs, forcing the stellar radius to in- 
flate in order to conserve flux dGough & T avler 1966). Recent 
attempts at modeling these effects have indicated that an en- 
hanced magnetic field is a plausible explanation, although the 
primary physical mechanism affecting the structure of t he star 
is still deba ted (Mullan & MacDonald 2001; Chabri er et all 
'2007';'MacDonald & Muflan 2012). 

Regardless of the precise physical mechanism, magnetic 
fields should betray their presence through the generation of 
magnetic activity in the stellar atmosphere. If magnetism is 
responsible for the observed inflated stellar radii, then cor- 
relations should be expected between individual stellar ra- 
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dius deviations and magnetic activity indicators (i.e., chro- 
mospheric Ha and Call H & K emission, coronal x-ray 
emission, etc.). Tantalizing evidence of such correlations 
has been reported p reviously by Locez-Morales (2007) and 
[Morales et al] JIOOl . 

However, recent evidence appears to stand in contrast 
with the current t heory. Two syst ems, LSPM Jl 112 +7626 
jlrwin etaljr201 1') and Kepler- 16 (iDovle et al.l 120111) . were 
discovered that have wide orbits with approximately forty- 
one day periods. Despite this, both appear to display discrep- 
ancies with stellar evolution models. In these systems, the 
component stars should be evolving individually with mutual 
tidal interactions playing a negligible role in the overall angu- 
lar momentum evolution. The stars should be spi nning down 
over time due to magnetic breaking processes ( jSkumanichI 
meaning the stars should not be as magnetically active 
compared with short period binary systems. The contrast is 
particularly evident for LSPM J 11 12+7626, where a rotation 
period of sixty-five days was detected via starspot modulation 
in the out-of-eclipse light curve. Gyrochronology suggests 
that t he system has an age of approximately 9 Gyr (Barnes 
12010') and implies further that the secondary is likely slowly 
rotating and should, therefore, not shown signs of strong mag- 
netic activity or an inflated radius. 

A third syst em also appears t o defy the current hypothe- 
sis. KOI- 126 (iCarter et alJl2"oTlh is a hierarchical triple sys- 
tem with two low-mass, fully convective stars in orbit around 
a 1.35 Mq primary. The two low-mass stars are orbiting each 
other with a period of 1.77 days. Therefore, they should 
show signs of inflated radii due to enhanced magnetic activ- 
ity. However, it has been shown that the two low-mass, fully 
convective stars were in agreement with model predictions 
when considering their su per-solar metallicity and the age of 
the higher mass primary (Feiden et al. 201 1). This_agreement 
was further confirmed by Spada & Demarque (201^. 

Metallicity has been proposed previously as a solution to 
the observed M R discrepanc ies, but for the case of single 
field stars (Berg er et alj|2006l). This was contradicted shortly 
thereafter by Lopez -MoralesI ( l2007h . most notably for DEBs. 
Although, we must consider that the radius discrepancy- 
metallicity correlation is severely complicated by the fact that 
metallicities of M dwarfs are notoriously difficult to deter- 
mine observationally. 

Finally, developments in light curve modeling of spotted 
stars has generated interesting results. The presence of large 
polar spots may alter the light curve analysis of DEBs by mod- 
ifying the eclipse profile. These modifications lead to 2 - 4% 



uncertainties in the derived stellar radii ( 


Morales et al.ll2010t 


iWindmiller et alJl20iaiKraus et al.ll201 1 


). Thus far, only two 



DEBs (GU Boo and CM Dra) have been thoroughly tested for 
their sensitivity to spots. Systematic uncertainties may there- 
fore dominate the error budget, casting a shadow of doubt on 
the observed radius discrepancies, which are often made ap- 
parent due to the minuscule random uncertainties. 

The uncertainties and developments outlined above have 
motivated us to reevaluate the current state of the low-mass 
MR relationship. In what is to follow, we use a large grid 
of theoretical stellar evolution isochrones in an effort to com- 
pare the low-mass models of the Dartmouth Stellar Evolu- 
tion Program (DSEP) with DEB systems that have well con- 
strained masses and radii. We then explore how potentially 
unaccounted for systematic uncertainties have the ability to 
create the appearance of discrepancies when neglected and 
mask real ones when considered. Section |2] will present the 



DEB sample followed by a description of the stellar models 
in Section[3] The isochrone grid and fitting procedures will be 
explained in Section |4] Results will be presented in Section 
|5]followed by a discussion of the implications of our findings 
in Section|6] We conclude with a brief summary of the entire 
study in Section|7] 

2. DATA 

Our selection criteria mimic those of [Torres et alj (120101 
hereafter TAG 10) in so far as we require the random uncer- 
tainties in the mass and radius measurements to be less than 
3%. No data was disqualified due to perceived data quality 
issues or the original author's attempts to constrain possible 
systematic uncertainties. We also applied the criterion that 
the DEB system must include at least one component with a 
mass less than 0.8 Mq. This cut-off in mass is used to desig- 
nate "low-mass" stars. 

The mass cut-off was selected for two main reasons. First, 
the affects of age and metallicity on the structure of low-mass 
main sequence stars are suppressed compared to stars with 
masses of approximately 1.0 M© or above. Overall, this al- 
lows for less flexibility in fitting models to the observations, 
providing a more critical analysis of the stellar evolution mod- 
els. Second, some of the largest discrepancies between obser- 
vations and models are seen in the low-mass regime. This is 
likely a consequence of the former reason: true discrepancies 
become more apparent as the models become less sensitive to 
the input parameters. 

Stars used in this study are listed in Table[T]along with their 
observationally determined properties and original references. 
Our final sample consisted of eighteen DEB systems for a to- 
tal of thirty-six stars. Six of these systems are taken from 
TAG 10 who reanalyzed the available data using a common 
set of reduction and parameter extraction routines in an effort 
to standardize the process. While the original references are 
cited, the parameters listed are those derived by TAG 10 whose 
results were similar to the original values. 

A majority of the systems, ten in total, were published after 
the release of the TAG 10 review. Six are drawn from the study 
performed bv lKraus et al.l (1201 Ih . three are products of recent 
results from the Kepler Sp ace Telescope mission ( Carter et aP 
1201 U iDoyleetal .1120111 Bass et al. in prep.), and the fi- 
nal post- TAG 1() syste m was discovered by the M-Earth sur- 
vey dlrwin et al.ll201 Ih . The remaining t wo systems from our 
samp l e were announced before TAG 10 dLopez-Morales et al] 
'2006^, ^pez-Morales & Shawll2007l) . however, they were not 
included in the review for reasons related to either data avail- 
ability or data quality. 

On e final note: the Kepler systems KO I- 126 (ICarter et alJ 
1201 Ih and Kepler- 16 (iDovIe etal.ll2011h were not analyzed 
in a similar manner to the rest of the double-lined DEB pop- 
ulation. They are DEB systems whose parameters were de- 
rived from Kepler photometry using a dynamical-photometric 
model (see Supporting Online Material from Carter et al] 
1201 ll) . However, they still satisfy the criteria for comparison 
with stellar evolution models and have been included for this 
reason. 

3. MODELS 

Models utilized in this study were computed with the 
Dartmouth Stellar Evolution Program (DSE Pfl a descen- 
dant of the Yale Rotating Evolution Code (iGuenther et alJ 

^ http://stellar.dartmouth.edu/~models/ 
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Table 1 

DEB systems with at least one low-mass component with precise masses and radii. 



Star 


Porb 


Mass 


Radius 


^rot 


Source ' 


Name 


(day) 


Wj) 




(day) 




UV Psc A 


0.86 


0.9829 ± 0.0077 


1.110 ± 0.023 




1 


UV Psc B 




0.76440 ± 0.00450 


0.8350 ± 0.0180 


0.80 




IM Vir A 


1.309 


0.981 ± 0.012 


1.061 ± 0.016 




2 


IM Vir B 




0.6644 ± 0.0048 


0.6810 ± 0.013 


1.31 




KID 6131659 A 


17.528 


0.924 ± 0.008 


0.8807 ± 0.0017 




3 


KID 6131659 B 




0.683 ± 0.005 


0.6392 ± 0.0013 






RX J0239. 1-1028 A 


2.072 


0.7300 ± 0.0090 


0.7410 ± 0.0040 




4 


RX J0239. 1-1028 B 




0.6930 ± 0.0060 


0.7030 ± 0.0020 






Kepler- 16 A 


41.08 


0.6897 ± 0.0034 


0.6489 ± 0.0013 




5 


Kepler- 16 B 




0.20255 ± 0.0007 


0.22623± 0.0005 






GU Boo A 


0.49 


0.61010 ± 0.00640 


0.6270 ± 0.0160 


0.49 


6 


GU Boo B 




0.59950 ± 0.00640 


0.6240 ± 0.0160 


0.54 




YY Gem A 


0.81 


0.59920 ± 0.00470 


0.6194 ± 0.0057 


0.87 


7 


YY Gem B 




0.59920 ± 0.00470 


0.6194 ± 0.0057 


0.82 




MG 1-506664 A 


1.55 


0.584 ± 0.002 


0.560 ± 0.005 




8 


MG 1-506664 B 




0.544 ± 0.002 


0.513 ± 0.009 






MGl-1 16309 A 


0.827 


0.567 ± 0.002 


0.552 ± 0.004 




8 


MGl-1 16309 B 




0.532 ± 0.002 


0.532 ± 0.004 






MGl-1819499 A 


0.630 


0.557 ± 0.001 


0.569 ± 0.002 




8 


MGl-1819499 B 




0.535 ± 0.001 


0.500 ± 0.003 






NSVS 01031772 A 


0.368 


0.5428 ± 0.0027 


0.5260 ± 0.0028 




9 


NSVS 01031772 B 




0.4982 ± 0.0025 


0.5088 ± 0.0030 






MGl-78457 A 


1.586 


0.527 ± 0.002 


0.505 ± 0.008 




8 


MG1-78457B 




0.491 ±0.001 


0.471 ±0.009 






MGl-646680 A 


1.64 


0.499 ± 0.002 


0.457 ± 0.010 




8 


MG1 -646680 R 

iVlvJ i U^UUOW 13 




0.443 ± 0.002 








MG 1-20563 16 A 


1.72 


0.469 ± 0.002 


0.441 ± 0.004 




8 


MG1-2056316B 




0.382 ± 0.001 


0.374 ± 0.004 






CU Cnc A 


2.77 


0.43490 ± 0.00120 


0.4323 ± 0.0055 




10 


CU Cnc B 




0.39922 ± 0.00089 


0.3916 ±0.0094 






LSPM 11 112-1-7626 A 


41.03 


0.3946 ± 0.0023 


0.3860 ± 0.0052 




11 


LSPM 111 12-1-7626 B 




0.2745 ± 0.0012 


0.2978 ± 0.0046 






KOI-126B 


1.77 


0.2413 ±0.0030 


0.2543 ± 0.0014 




12 


KOI-126C 




0.2127 ± 0.0026 


0.2318 ±0.0013 






CM Dra A 


1.27 


0.23102 ±0.00089 


0.2534 ± 0.0019 




13 


CM Dra B 




0.21409 ± 0.00083 


0.2398 ±0.0018 






(1) Popper' ('19971; 


(2) Morales et al. 1 2009b); (3) 


G. Bass et al. 


(in prep.); (4) 



Lopez -Morales & Shaw (2007); (5) Doyle et al. (2011); (6) Lopez-Morales & Ribas (2005); 



(7)ToiTes & Ribas (2002); (8) Kraus et al. (2011 J; (9) iLope z-Morales et al. (2006); (10)|Ribai 
(2003); (11) Irwin et al. (201 1); (12) Carter et al. ROTll) ; (13) Morales et al. (2009a); 

Il992h . Physics included in the mod els have been 

described extensively in the literature (IChab over et al.l 

2001; Biork & ChaboyiB |200i iDotter et alj 12007. 2008; 
Feiden et al.i201 1 ). Below, we will summarize the physics in 



corporated in the code that are pertinent for the present work. 

Of particular importance for work on low-mass stars is the 
equation of state (EOS). In general, the EOS in DSEP is mass 
dependent. Stars with masses above O.SM© are well repre- 
sented by an ideal gas EOS with an appropriate Debye-Hiickel 
correct ion to account for ion charge shielding by free elec- 
trons (Chabover & Kim 1995). Below O.8M0, non-ideal con- 
tributions to the EOS become non-negligible and must be con- 
sidered. In this case, we use the FreeEOSQ in the EOS4 con- 
figuration. FreeEOS allows for the treatment of an arbitrary 
heavy element abundance, providing DSEP with the flexibil- 
ity to more reliably calculate models with super-solar abun- 
dances. 

Conditions in the outer, optically thin layers of low-mass 
stars precludes the use of grey atmosphere approximations 
and t he use of radiative T(t) relations ( Chabri er & Baraff^ 
120001 and references therein). As such, it is critical to apply 
the results of non-grey model atmospheres to the definition 
of surface boundary conditions. For consistency, calculations 

' by Alan Irwin, p5p://freeeos .sourceforge .net | 



across all mass regimes in this study define the surface bound- 
ary conditions using t he PHOENIX AMES -COND model at- 
mospheres (Hauschil dt et aljri999allbh . The atmospheres are 
attached at T = Teff by interpolating in tables generated by 
PHOENIX. Interpolation is performed in two variables, \ogg 
and Pgas, in order to define the temperature at the surface of the 
star. By attaching the atmospheres atT = T^ff, we eliminate 
the need to select an exact co nvective mixing-len gth value in 
the atmosphere calculations (iBaraffe et al.lll997b . However, 
this only applies to stars above O.2M0. Below that threshold, 
it becomes imperative to match the surface boundary condi- 
tions at a deeper optical depth in order to maintain an adia- 
batic atmosphere profile (Ch abrier & Baraff5l2000t) . We do 
not concern ourselves with this, for now, as none of the stars 
in our sample have a mass below the threshold. 

Stars above about O.SSM© begin to develop radiative cores 
in our model set, with the size of the surface convection zone 
shrinking as mass increases. With this in mind, accurate ra- 
diative opacities ar e required. We use the O PAL high tem- 
perature opacities dlglesias & Rogers! 1 1996h in conjunctio n 
with the low temperature opacities of 'Fergu son et al.l (l2005h . 
The Ferguson low temperature opacities are also utilized in 
the PHOENIX model atmospheres, providing consistency be- 
tween the atmosphere calculations and the stellar envelope in- 
tegration within DSEP. 
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As stars develop radiative cores it becomes increasingly 
important to include diffusion ph ysics. Helium an d heavy 
element diffusion are included ( Thoul et alj 11994 !). unless 
the star is fully convective. Fully convective stars are as- 
sumed to be completely and homogeneously mixed through- 
out since the convective timescale dominates the diffusion 
timescale. When diffusion is calculated, a turbulent mix- 
in g term is included in the diffusion equations as prescribed 
by iRichard et alj (l2005h . The adopted reference temperature 
characterizing the efficiency of turbulent mixing is log(rref) = 
6.0 following the analysis of Korn et al. (2007). 

Most of the stars in our sample may be considered rapid ro- 
tators. Therefore, one naturally wonders whether we should 
consider the effects of rotational deformation in our mod- 
els. As a first approximation, we applied Chandrasekhar's 
analysis of slowly ro tating polytropes to our DEB sample. 
IChandrasekhar' ("1933) derived an expression for the stellar 
oblateness analytically. He defined slowly rotating to be when 



27iGp, 



< 1 



(1) 



where co is the stellar angular velocity and is the central 
mass density. Assuming the stars in our sample are spin-orbit 
coupled, we find that all but two stars have % < lO^'*. The two 
exceptions (one DEB system) have % w 10^^; a consequence 
of their short orbital period of 0.368 days. Hence, all of the 
stars in our sample satisfy the slowly rotating approximation. 

Based on the assumed polytropic index, n, Chandrasekhar 
derived that the deviation from sphericity (oblateness) could 
be approximated as 



' eq 



5.79% 
9.82X 
41.81X 



for n = 1.5 
for n = 2.0 
for n = 3.0 



with !f being the oblateness, and r^q and rpoie are the equato- 
rial and polar radius, respectively. We approximated the poly- 
tropic index for each star using the interior density profile pre- 
dicted by DSEP. The density profile for low-mass stars only 
slightly deviates from the polytrope prediction over the inner 
98% of the star (by radius). Specifically, we found that be- 
low 0.4 Mq, the models were best represented by an n = 1.5 
polytrope. Above 0.4 Mq but below 0.65 Mq, an n = 2.0 poly- 
trope was appropriate. For all masses greater than 0.65 Mq, 
we defaulted to assuming an n = 3.0 polytropic index. 

For all of the systems in question we found < 0.001, 
except for the one system that had the higher % value. This 
system, NSVS 01031772, is slightly deviating from spheric- 
ity with ^ « 0.011. We feeljustified in neglecting the physics 
of rotation in our stellar models, with the caveat that attempt- 
ing to probe model precisions below 1% for NSVS 01031772 
will likely require a more detailed treatment of rotational de- 
formation. 

Finally, we introduced the capability for DSEP to compute 
the characteristic convective overt urn time. Our implemen - 
tation is similar to the method of Kim & Demarqud (119961) , 
who calculate the "local" convective overturn time at a dis- 
tance above the base of the convection zone equal to one 
half of the mixing-length. This particular location was cho- 
sen based on the assumption that the tachocline (radiative- 
convective zone in terface) is the source region of the stellar 
magnetic dynamo (lParkerlll975h . However, using the stellar 
tachocline as our magnetic field source location for a fully 
convective star would be nonsensical. Therefore, in the case 



Table 2 

Solar Calibration Parameters 



Parameter 



Adopted 



Model 



Age (Gyr) 

Me(g) 
Rq (cm) 
L@ (erg s"') 

{Z/X)s,irf 



4.57 
1.9891 xlO'3 
6.9598 xlO'" 
3.8418 xlO^^ 
0.713±0.001 
0.0231 



log(R/Re) = 8 X 10-^ 
log(L/L;:0 = 2xlO-'* 

0.714 

0.0230 



th at the star i s conv ective throughout, we utilize the results 
of iBrowningI (l2008h as a first approximation to the magnetic 
field source location. Browning found that the magnetic field 
strength within a fully convective star was at a maximum at 
a depth of 85% of the stellar radius. In accordance with this 
result, we compute the convective overturn time at one half 
the mixing-length above this location. 

4. ISOCHRONE FITTING 

4.1. Solar Calibration 

Carrying out a detailed comparison of theoretical models 
with observational data first requires the determination of a 
precise solar calibration configuration. This calibration is per- 
formed in order to determine the appropriate proto-solar he- 
lium mass fraction (Y), heavy element mass fraction (Z), and 
the convective mixing-length {<X mlt)- Using the he a vy ele - 
ment abundance composition of iGrevesse & Sauvall (Il998l) . 
a X.OMq model w as evolved to the solar age (4.57 Gyr, 
iBahcall et al.l 120051) at which point the model was required 
to reproduce the solar radius, solar luminosity, observed ra- 
dius to the base of the convection zone, and the solar photo- 
spheric (Z/X). The final set of parameters required to match 
the solar properties were F,>„-, ~ 0.27491, Z,,,,-, ~ 0.01884, and 
CMLr = 1.938. Corresponding solar model parameters are 
listed in Table |2] 

4.2. Isochrone Grid 

Isochrones were computed for a wide range of age and 
metallicity values. The parameter space was defined to en- 
compass a vast majority of stars typical of the local galac- 
tic neighborhood. Seven ages and seven scaled-solar metal- 
licities were adopted for a total of forty-nine individual 
isochrones. The sets of values used in this study were: [Fe/H] 
= {-1.0, -0.5, -0.3, -0.1, 0.0, H-0.1, H-0.2} dex and age = {0.3, 
1.0, 3.0, 5.0, 6.0, 7.0, 8.0} Gyr. For completeness, we also 
generated a set of isochrones that employed a smoothly vary- 
ing convective mixing-length. However, for clarity, we defer 
the discussion of these models and relegate the information to 
AppendixlAl 

4.3. Fitting Procedure 

Judging agreement between observations and models was 
performed on a system-by-system basis, as opposed to fit- 
ting individual stars. Critical to the process was ensuring 
that both components of a given system were consistent with 
isochrones of a common age and metallicity. For each ob- 
ject, corresponding model radii were derived by linearly in- 
terpolating in each isochrone using the observationally de- 
termined mass. A linear interpolation scheme was sufficient 
since the mass resolution along the isochrones was small 
(AM = O.O2M0). 

Relative errors between the model radii and the observa- 
tionally determined radius were then calculated. The relative 
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Table 3 

Age and [Fe/H] priors 



DEB System 


Age (Gyr) 


[Fe/H] 


Status 


UVPsc 


7.9 




Rejected 


IMVir 




-0.3±0.3 


Rejected 


Kepler- 16 




< 0.0 


Accepted 


GU Boo 


< 1.0 




Rejected 


YYGem 


0.4 ±0.1 


+0.1±0.2 


Accepted 


CU Cnc 


0.3 ±0.1 




Rejected 


KOI-126 


4.0 ± 1.0 


+0.15±0.8 


Accepted 


CMDra 


4.0 ± 1.0 


<0.0 


Accepted 



error was defined as 



dR 

Robs 



obs 



^ mod el 



Robs 



(2) 



and will be presented as such throughout the rest of the paper. 
Formal agreement was determined by analyzing the number 
of standard deviations outside of the accepted range our model 
radii were located. Explicitly, 



#Or = 



Robs Rjnodel 
Or 



(3) 



where Or is the random uncertainty of the observational ra- 
dius. 

Different levels of compatibility were assigned to each sys- 
tem based on the individual agreement of each component 
in the system. For instance, if an isochrone matched both 
components at a common age within the la limits set by the 
quoted random uncertainties, the system as a whole was des- 
ignated as being "fit" by the models. Only within the la level 
was a system considered to be accurately represented by the 
models. The analysis was applied to two sets of data, which 
differ in the adopted radius uncertainties (formal quoted un- 
certainties or fixed 3% uncertainties) that will be justified 
later. 

The process detailed above generated a list indicating the 
level of agreement between each isochrone and each DEB 
system. Narrowing the list to a single "best fit" isochrone 
involved minimizing the root mean square deviation (RMSD) 



RMSD 




(4) 



where the sum is performed over the components of the DEB 
system. 

4.4. Age & Metallicity Priors 

DEBs that have been well studied have additional con- 
straints that allow us to restrict the set of isochrones used in 
the fitting procedure. Specifically, our sample contains eight 
systems that have the added constraint of either an estimated 
age or metallicity, and in a couple cases, both. Before ac- 
cepting the quoted age and metallicity priors, however, we 
performed a qualitative analysis to ensure that the estimates 
we adopted were reliable. A summary of our analysis is pre- 
sented in Appendix |B] We judged four of the eight systems 
with quoted age or metallicity priors to have been determined 
reliably. Adopted priors are listed in Table |3] 

5. RESULTS 
5.1. Direct Comparison 
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Figure 1. Relative en'ors between stellar evolution models and observation- 
ally determined radii. Reliable age and metallicity priors were accounted 
for in the statistical analysis. (Top) The adopted measurement uncertainties 
are the cited observational random uncertainties derived from the light curve 
fitting procedure. (Bottom) A fixed 3% uncertainty is adopted to represent 
possible systematic uncertainties (i.e., starspots). 

Figure [T] demonstrates that, in general, our models reduce 
the observed radius discrepancies to below 4% for 92% of the 
stars in our sample. Only a few outliers are seen to be largely 
discrepant. Across the entire sample, we find a mean absolute 
error of 2.3%. This broadly represents a factor of two reduc- 
tion in the previously cited radius discrepancies. The age and 
metallicity of the best fit isochrone for each DEB system is 
listed in Table |4] along with the level of agreement between 
the best fit isochrone and the individual stars comprising each 
DEB system. 

Contrasting with previous studies, we do not observe an 
overwhelming systematic trend of the models grossly under 
predicting stellar radii. Instead, systems discrepant by more 
than 5% represent an exception to the broad agreement be- 
tween the models and observations. While the agreement 
is admittedly not perfect, it is apparent that most observed 
discrepancies must now be judged according to the precision 
with which their radii were measured. As we will see in the 
next subsection, systematic uncertainties have the potential to 
blur our interpretation of the agreement between models and 
observations, simply due to the factor of two improvement 
described above. 

Ignoring systematic uncertainties, for now, we are still con- 
fronted with the need to explain both the slightly discrepant as 
well as the largely discrepant radii. Discussion and comments 
pertaining to plausible explanations of the inflated radii will 
be addressed later in Section|6] 

5.2. Fixed Uncertainties 

Recent work in modeling eclipse profiles of DEB sys- 
tems has led to the realization that large, polar starspots 
can potentially affect observationally derived radii and sub- 
sequently domin ate over the quoted ra ndom uncerta inties. 
IWindmiller et al.' (GU Boo; ,20 IQ) and [Morales et all (CM 
Dra: 120 lol) have found that systematic uncertainties may be 
present in DEB radii on the order of 2 - 4%. This is to be 
compared with quoted random uncertainties that are typically 
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Table 4 

Best fit isochrone with a solar caHbrated mixing-length. 



Quoted Fixed 3% 



Star Name 


Age 


[Fe/H] 


Sff/Robs 




Fit 


Age 


[Fe/H] 


oR/Robs 




Fit 


UV Psc A 


8.0 


-0.10 


-0.0285 


-1.375 


No 


8.0 


-0.10 


-0.0285 


-0.950 


No 


UV Psc B 






0.1031 


4.785 








0.1032 


3.439 




IM Vir A 


7.0 


-0.10 


-0.0037 


-0.248 


No 


7.0 


-0.10 


-0.0038 


-0.125 


No 


IM Vir B 






0.0408 


2.136 








0.0408 


1.359 




KID 6131659 A 


3.0 


-0.50 


-0.0014 


-0.750 


Yes 


3.0 


-0.50 


-0.0014 


-0.047 


Yes 


KID 6131659 B 






0.0004 


0.200 








0.0004 


0.013 




RX J0239. 1-1028 A 


8.0 


-0.10 


0.0297 


5.510 


No 


8.0 


-0.10 


0.0297 


0.991 


Yes 


RX J0239. 1-1028 B 






0.0285 


10.007 








0.0285 


0.949 




Kepler- 16 A 


1.0 


-0.10 


0.0004 


0.210 


No 


0.3 


-0.10 


0.0099 


0.330 


Yes 


Kepler- 16 B 






0.0299 


5.510 








0.0225 


0.749 




GU Boo A 


8.0 


-0.10 


0.0292 


1.144 


No 


8.0 


-0.10 


0.0292 


0.973 


No 


GU Boo B 






0.0414 


1.616 








0.0414 


1.381 




YY Gem A 


0.3 


-0.10 


0.0813 


8.834 


No 


0.3 


-0.10 


0.0813 


2.710 


No 


YY Gem B 






0.0813 


8.834 








0.0813 


2.710 




MG 1-5 06664 A 


1.0 


-0.10 


0.0035 


1.951 


No 


1.0 


-0.10 


0.0035 


0.116 


Yes 


MG 1-5 06664 B 






-0.0031 


-1.612 








-0.0031 


-0.105 




MGl-1 16309 A 


7.0 


-0.50 


-0.0084 


-1.162 


No 


7.0 


-0.50 


-0.0084 


-0.281 


Yes 


MGl-1 16309 B 






0.0150 


1.993 








0.0150 


0.500 




MGl-1819499 A 


8.0 


-0.50 


0.0296 


8.425 


No 


7.0 


-0.10 


0.0404 


1.348 


No 


MGl-1819499 B 






-0.0528 


-8.792 








-0.0408 


-1.360 




NSVS 01031772 A 


8.0 


-1.00 


-0.0021 


-0.396 


No 


8.0 


-0.50 


-0.0134 


-0.447 


No 


NSVS 01031772 B 






0.0412 


6.991 








0.0395 


1.317 




MGl-78457 A 


5.0 


+0.20 


-0.0008 


-0.049 


Yes 


5.0 


+0.20 


-0.0008 


-0.026 


Yes 


MGl-78457 B 






-0.0001 


-0.006 








-0.0001 


-0.004 




MG 1-646680 A 


1.0 


+0.20 


-0.0226 


-1.725 


No 


1.0 


+0.20 


-0.0226 


-0.755 


Yes 


MG 1-646680 B 






0.0184 


1.309 








0.0184 


0.613 




MGl-2056316 A 


3.0 


+0.20 


-0.0074 


-1.632 


No 


3.0 


+0.20 


-0.0074 


-0.247 


Yes 


MG1-2056316B 






0.0067 


1.248 








0.0067 


0.222 




CU Cnc A 


8.0 


+0.20 


0.0233 


1.832 


No 


8.0 


+0.20 


0.0233 


0.777 


Yes 


CU Cnc B 






0.0020 


0.084 








0.0020 


0.067 




LSPMJ1112-H7626 A 


8.0 


+0.20 


-0.0005 


-0.035 


No 


8.0 


+0.20 


-0.0005 


-0.016 


Yes 


LSPM 111 12+7626 B 






0.0290 


1.875 








0.0290 


0.966 




KOI- 126 B 


1.0 


+0.10 


-0.0011 


-0.205 


Yes 


3.0 


+0.20 


-0.0036 


-0.120 


Yes 


KOI-126C 






-0.0003 


-0.053 








-0.0004 


-0.013 




CM Dra A 


5.0 


0.00 


0.0316 


4.221 


No 


5.0 


0.00 


0.0316 


1.055 


No 


CM Dra B 






0.0360 


4.798 








0.0360 


1.200 





on the order of 1% or les s. 

Similarly, iKraus et aTl (1201 Ih included an estimate of sys- 
tematic uncertainties for their radius measurements and found 
typical uncertainties on the order of 2 - 3%. With the ra- 
dius uncertainties potentially dominated by often unquoted 
systematics, we were curious to see what effect such uncer- 
tainties would have on comparisons with theoretical models. 
Quite obviously, observed discrepancies would be reduced by 
the level of systematic uncertainty. However, it has not been 
clear whether those systematics have the ability to completely 
mask the observed radius residuals, particularly when com- 
bined with isochrones that cover a wide age-metallicity pa- 
rameter space. Since the radius uncertainty cut-off for our 
DEB sample was set at 3%, we elected to adopt a fixed 3% 
radius uncertainty in our analysis to mimic potential system- 
atics. 

Demonstrated in Figure[T]B is the effect of including a fixed 
3% radius uncertainty. The overall distribution of points ap- 
pears to be similar to the previous case, except that the larger 
uncertainty enables more systems (eleven in total) to be con- 
sidered fit by our analysis (see Table |4]l. This result is clearly 
expected when introducing larger error bars, as mentioned 
above. 

What is important to realize is that we are now presented 
with a scenario where the systematic uncertainties mask our 
ability to draw firm conclusions about whether the observed 
radius discrepancies are inherently real or merely a conse- 
quence of neglected uncertainties. Such ambiguity was pre- 



viously not present since the observed radius residuals were 
substantially larger than any potential systematic uncertain- 
ties. 

5.3. Peculiar Systems 

Before continuing with a further interpretation of our re- 
sults, we wanted to briefly comment on several individual 
DEB systems. These systems stood out in our mind as worth 
remarking upon separate from the ensemble. 

Our models never fit the UV Psc system to a coeval age. 
While we were able to fit UV Psc A, the secondary compo- 
nent was always found to have an 8 - 10% larger radius than 
wou ld be expected from stellar evolution theory, consistent 
with lPoppe^ (11997b . Both stars may realistically be discrepant 
with the models, however, if the age is truly younger than 
8 Gyr. Constraining the actual magnitude of the discrepan- 
cies is difficult without observational age and metallicity esti- 
mates. Fittingly, the system is known to be very active based 
on spectroscopic a nalysis of Ha cores, m ulti-epoch photo- 
metric monitoring dKiurkchieva et al.''2005 '), and the derived 
x-ray luminosity (this work, Section l6.2.2b . Further observa- 
tions of U V Psc to constrain its metallicity and to provide data 
for a more rigorous starspot analysis would be a worthwhile 
endeavor 

MGl-1 89499, characterized by iKraus et alj ( 120111) . was 
also never found to be in total agreement. The primary and 
secondary were found to deviate by 3.0% and -5.3%, respec- 
tively. This appeared to be a cause for concern, but Kraus et 
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al. were one of the few groups to provide an estimate of po- 
tential systematic uncertainties. Uncertainties in the primary 
star's radius were elevated to 4.6% with a comparatively large 
uncertainty of 3.4% in the secondary's radius. Applying these 
uncertainties and rerunning the isochrone fitting procedure al- 
lows the predictions of the models to nearly fall inline with 
the observations. An age of 8 Gyr was again derived but with 
a new metallicity of [Fe/H] = -0.3. This isochrone yielded 
relative radius errors of 4.8% and -2.9% for the primary and 
secondary, respectively. Hence, the secondary is now consid- 
ered fit by the models and the primary is only 0.2% outside 
the bounds of uncertainty. 

The last two systems we would like to discuss are two 
that have long pose d problems for model ers: YY Gem and 
CM Dra. YY Gem ( iTorres & Ribasll2002l) is an equal mass, 
equal radius DEB system that is effectively represented by a 
single point at (0.6, 0.08) in Figure [T] Beginning with the 
low-mass models of Hoxie ( 1970), YY Gem has never been 
adequately reproduced by stellar models. 

The difficulty with YY Gem seems to be in the estimated 
age of 400 Myr, determined by its association with the Cas- 
tor AB quadruple system. An older age is logically preferred 
by the models considering the notably inflated radius. Inter- 
estingly, this would have led previous studies to underesti- 
mate the observed radius discrepancy. Assuming either a 1 or 
5 Gyr isochrone naturally results in models with larger radii 
than those computed with an age of 400 Myr. Further study 
of this system, particularly its distance (Section r6.2.2b and as- 
sociation with Castor AB (see Appendix|B]i, will be beneficial 
toward fully understanding the nature of its inflated radii. 

Finally, CM Draconis is another constant thorn in theo- 
reticians sides. Here, the primary and secondary deviate 
by 3.2% and 3.6%, respectively. CM Dra is one of the 
most studied DEB systems and has very well constrained 
physical properties measured from data spanning multiple 
decades ([Morales et al. 2009a). Our models indicate that a 
solar composition is preferrecfl However, the composition of 
CM Dra, while not known positively, is very likely subsolar 
dViti et al.ll 19971 IIOO ?, 'Morales et al."2009a';^Kuz "netsov et all 
I2012^ ■ meaning our models are probably more discrepant than 
indicated by this study. Moral es et al.l (12010) found that large 
polar spots produce systematic uncertainties in the radius 
measurements of around 3%, potentially bringing CM Dra 
more in line or more out of line with model predictions. For 
the time-being, this system will continue to test our knowl- 
edge of stellar evolution. 

6. DISCUSSION 

Data presented in this study hint at two competing expla- 
nations for the occurrence of the differing model and ob- 
servational MR relations. It is entirely plausible that stel- 
lar evolution models are not incorporating key physical pro- 
cesses that can account for the observed discrepancies. This 
view is not new and has always accompanied discussions 
of low-mass models (Hoxie 1970, 1973; ChabrierABaraffd 
ll997HBaraffe et al.ll997L l998). Typically cited is our incom- 
plete knowledge dealing with the complex array of molecules 
present in M-dwarf atmospheres as well as the lack of struc- 
tural changes induced by a large-scale magnetic field. In- 
cluded in the latter are both the effects on convective energy 

^ A super-solar metall icity is actually favored, as was also noted by 
ISpada & Demarau3 j2012D . but the application of the metallicity prior in our 
analysis restricted the models to only solar or subsolar compositions. 



transport and the emergence of spots on the stellar photo- 
sphere. 

The other scenario is one in which the neglect of systematic 
uncertainties is driving the apparent discrepancies. We have 
shown that by allowing for realistic variation in age and metal- 
licity of the models, the radius residuals are of the same mag- 
nitude as the potential systematic uncertainties. Note, this is 
without any modification to the solar calibrated mixing-length 
or the inclusion of any non-standard physics. Previous stud- 
ies have indicated that systematics may help alleviate the size 
of the radius residuals, but have required additional modifica- 
tions to the models in order to fully reconcile models with ob- 
servations. It is now clear that we must work to constrain and 
minimize systematic uncertainties in observations of DEBs to 
allow for them to provide an accurate test of stellar evolution 
models. 

6.1. Radius Deviations 

The MAE between our models and the observations was 
2.3%, a factor of two improvement over the canonical mini- 
mum of 5%. We found deviations of no more than 4%, with 
the exception of a few stars, instead of ubiquitous 5-15% 
errors, as is often quoted. Despite improving the situation, 
panel A of Figure [TJillustrates that the models are still unable 
to fully reproduce the observed stellar radii. 

Accepting the factor of two improvement presented in Sec- 
tion l5.1l the paradigm of broad disagreement between models 
and observations is shifted to one where agreement is broad, 
and large discrepancies are an exception. With the radius de- 
viations typically less than 4%, an evaluation of the system- 
atic errors becomes imperative. Formerly, systematic uncer- 
tainties of about 4% were incapable of relieving radius devia- 
tions greater than 5%. Stellar evolution models still appeared 
to disagree with DEB observations even after the inclusion of 
systematic errors. 

The reason for the factor of two improvement is twofold. 
First, we have calculated models with a finer grid of metallic- 
ities. DSEP utilizes an EOS that enables models to be more 
reliably calculated for super-solar metallicities, allowing for a 
greater range of stellar compositions to be considered. Low- 
mass stellar models with super-solar metallicity have previ- 
ously been unavailable for comparison with low-mass DEB 
data. Thus, the ability to extend our model set to super-solar 
metallicities allows for more flexibility in attempting to match 
the observed properties of DEB components. 

Second, a far larger number of low-mass DEBs with pre- 
cisely measured radii were available to us as opposed to pre- 
vious studies. Before the publication of TAGIO, there were 
only eight systems that met the criteria necessary to accurately 
constrain stellar evolution models. Following the TAGIO re- 
view, the total number of systems that met the necessary cri- 
teria more than doubled with the addition of ten newly char- 
acterized DEBs. These additional systems appear to be more 
in line with the results of standard stellar evolution theory. 
However, well-known discrepant systems remain noticeably 
discrepant and still require further explanation. 

We must now ask, "what belies the current discrepancies 
between models and observations?" Figure [T] favors the hy- 
pothesis that non-standard physics are absent from current 
stellar evolution models. Our larger data set allows us to no- 
tice that stars of similar masses from different DEB systems 
appear to be discrepant at varying levels, an effect a single 
set of standard models cannot correct. However, this is con- 
tingent upon the accuracy our age and metallicity predictions. 
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Figure 2. Observed coiTelation between radius deviations and DEB orbital peiiods. Maroon asterisks ai'e stai's with known x-ray flux measurements. Dashed 
vertical lines represent the various period cuts used in our statistical analysis to divide the sample into two subsamples. (Left) The full range of periods present in 
the current sample. (Right) Highlighting the short period regime as it encompasses a majority of the DEB systems. 



Until we have better empirical age and metallicity estimates 
for the various systems, it is too difficult to ascertain the true 
level of discrepancy for any individual star 

The efficiency of convective energy transport is of greatest 
interest. It is possible that convection is naturally inefficient. 
Although, we gather from Figure [T] that suppression of con- 
vective energy transport must be tied to a stellar property that 
is largely independent of mass. Simple parametrization of the 
suppression of convection is too uniform over a given mass 
regime to fully account for the observed differences in stel- 
lar radii for stars with similar masses (see also Appendix [All. 
Any effort, either theoretical or observational, to constrain the 
physics of convection in low-mass stars will lend crucial in- 
sight. 

The most favored option, is that convection is not intrin- 
sically inefficient, but that a la rge-scale magnetic field acts 
to suppress convective motions (iMulla n & MacDonaldll200ll: 
iChabrier et all 120071: iMacDon ald & Mullan| |2012|) Stellar 
evolution models self-consistently incorporating the effects of 
large-scale magnetic fields will help on this front. Observa- 
tions of cool-star magnetic field strengths and topologies will 
then provide a means of judging the validity of any new mod- 
els. 

One final hypothesis is that starspots may affect the struc- 
ture of stars and ge nerate the inflated radii we observe. 
IChabrier et al.l ( 120071) investigated such a possibility by arti- 
ficially reducing the total stellar bolometric luminosity in an 
effort to mimic the effects of spots. They found that radius 
discrepancies were relieved with their parametrization. How- 
ever, it is still not apparent whether spots reduce the total 
bolometric flux or if they locally shift flux to longer wave- 
lengths (Jackson et al. 2009), preserving the total luminosity. 
Ultimately, if starspots are required in stellar evolution mod- 
els, their inclusion is necessitated in the analysis of DEB light 
curves. 

Quantifying the effect starspots may have on observed DEB 
light curves is extremely difficult. Obtaining accurate knowl- 
edge of the total surface coverage of spots, the total number of 
spots, their individual sizes, temperature contrasts, and their 
overall distribution on the stellar surface is nearly an impos- 
sible task given only a light curve. From a theoretical per- 
spective, finding a proper parametrization to mimic the effect 
of spots on a three dimensional volume within the framework 



of a one dimensional model provides its own complications. 
Currently, the only feasible method to include spots, is to in- 
clude their potential effects on the radius measurement uncer- 
tainties. 

Unfortunately, while the inclusion of fixed 3% radius un- 
certainties in our analysis was able to alleviate many of the 
noted radius discrepancies, it also created a situation where 
the measurement uncertainty was on the order of the typical 
radius deviation. We are presented with a case where the ob- 
servations are no longer effective at testing the models and the 
manifestation of most radius discrepancies can be attributed 
to under estimated error bars. At this point, we require obser- 
vations that have been rigorously vetted for potential system- 
atic uncertainties and are able to still provide mass and radius 
measurements to better than 2%. 

6.2. Radius-Rotation-Activity Correlations 

Direct measurements of low-mass magnetic field strengths 
are rare, especially among fast rotating stars (Reiners 2012). 
Without a direct measure of the magnetic field strength, we 
are forced to rely on indirect measures to probe correlations 
between stellar magnetism and the appearance of inflated stel- 
lar radii. Ideally, these indirect measures are intimately con- 
nected with the dynamo mechanism, thought to generate and 
maintain stellar magnetic fields, or are the product of mag- 
netic processes in the stellar atmosphere. The preferred indi- 
rect measures are typically stellar rotation or the observation 
of magnetically driven emission (Ha, Call H and K, X-ray). 

6.2.1. Rotation 

Typically, low-mass DEBs are found in tight, short period 
orbits (< 3d; see Figure|2|. Tidal interactions spin-up the indi- 
vidual components and allow them to remain rapidly rotating 
throughout their life cycle. Stellar dynamo theory dictates that 
the large-scale magnetic field strength is tied to the rotational 
properties of a star (i.e., a rapidly rotating star should be more 
magnetically active than a comparable star that is slowly ro- 
tating; Parker 1979; Reiners et al. 2012), providing a natural 
starting point for our investigation. 

We performed two independent statistical tests on the dis- 
tribution of radius residuals as a function of the orbital pe- 
riod (Porb)- Our primary objective was to determine whether 
rapidly rotating systems produce, on average, larger radius 
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Figure 3. The theoretical Rossby number, Ro = Pjot/tconv, versus the relative 
radius error. Ro is tied directly to the theoretical stellar dynamo mechanism 
and is empirically related to the ratio of a star's x-ray to bolometric luminosity 
(a magnetic activity indicator). Asterisks in maroon are stars with known x- 
ray flux measurements. 

deviations than systems perceived to be slow rotators. Ro- 
tational periods were assumed to be synchronous with the 
orbital period unless a separate value for the rotational pe- 
riod was cited in the literature. The statistical tests performed 
were a Kirmogorov-Smirnov (K-S) test and another whereby 
we tested the probability of obtaining a given distribution of 
residuals via a Monte Carlo method. Both tests were per- 
formed on the observed difference in mean absolute error 
(MAEj3 between two data bins. The data bins were divided 
at preselected values of Porb^ identified visually as vertical 
dashed lines in Figure |2] 

Comparing the radius deviations with the rotational periods, 
we find no evidence of any dominant correlation. Figure |2] 
displays the residual data as a function of the orbital period, 
with frame A showing the full range of observed periods and 
frame B highlighting the "short period" regime. The correla- 
tion of radius deviations with orbital period has been studied 
previously (e.g., Kraus et al. 201 1) where a significant differ- 
ence between the two bins was observed around Porb = 1-5 
days. We performed the statistical tests using three values for 
the orbital period (l.Od, 1.5d, and 2.0d) that defined the two 
period bins. 

We confirm the results of lKraus et alj (1201 Ih and find a 3 . 1 a 
difference in the distribution of radius deviations around 1.5d. 
Systems with Porb < 1 -Sd had a MAE of 3.4% while the longer 
period systems had a MAE of 1.0%. While this is tantalizing, 
we can not necessarily attribute any physical significance to 
this particular division. We should expect this difference to 
be present for any two subsamples. However, we fail to find 
any evidence for a statistically meaningful difference when 
we divide the subsamples at l.Od and 2.0d. Therefore, the sig- 
nificant difference noted at 1 .5d is likely a spurious statistical 
resulfl Inclusion of more long-period DEBs will be instru- 
mental in providing a robust conclusion. Until those systems 

^ We selected the MAE over the RMSD as a measure of the mean radius 
deviation of a given ensemble in order to reduce the weight of any individual 
outlier in the final mean. 

^ Kra us et al] 420 llh posit the difference may actually be a by-product of 
the DEB light curve analysis methods. Providing a further examination is 
outside the scope of this study. 



are discovered, there does not appear to be a physically mean- 
ingful explanation for why the divide should be made at 1.5d, 
but then also not hold for a division at l.Od. 

The rotational (or orbital) period is not necessarily the most 
appropriate proxy for the magnetic field strength or poten- 
tial magnetic activity level that we could select. It would be 
ideal for the rotational parameter to have some connection to 
intrinsic stellar properties. For instance, rotational velocity 
normalizes the rotational period to the stellar radius, provid- 
ing a distinction between two main sequence stars of different 
masses that may have similar rotation periods. Optimally, the 
rotational variable in question would also provide a direct link 
to either observable magnetic activity or a theoretical descrip- 
tion of stellar magnetism. 

Accordingly, we advocate the use of the Rossby number 
(hereafter Ro). Ro is defined as the ratio of the rotational pe- 
riod of the star to its convective overturn time, Ro = Prot/tconv, 
and measures the strength of the Coriolis force acting on the 
vertical motion of convection cells. The dimensionless quan- 
tity Ro appears directly in standard mean-field dynamo theory 
(a-CO dynamo; lParkei|[r979in and is intimately related to the 
ratio of the stellar x -ray luminosity to bolometric luminosity 
('Wright et al."20lj). The latter quantity has been shown to be 
a strong indicator of a stellar corona heated to over 10^ K by 
magnetic activity (IVaiana et alJ[T98l IPailavicini et alJ[T98l 
iNoves et al.|[T984h . 

One cut in Ro was performed at Ro = 0. 1 , illustrtated in Fig- 
ure [3] The selection of Ro = 0. 1 approximately corresponds 
to Rosat, or the value of Ro associated with an observed sat- 
uration of the stellar dynamo apparent in the ratio of the stel- 
lar x-ray lumino sity to the bolometric luminosity (i.e., coro- 
nal saturation; Wrig ht et alJl201 If) . Intuitively, this suggests 
that all points with Ro values less than 0.1 are, presumably, 
sufficiently active so as to display inflated radii. Given our 
current understanding of coronal saturation, it is difficult to 
ascertain how strong of correlation is expected to exist. With 
that said, we assume that stars with very low Ro values should 
show at least a marginal degree of inflation compared to stars 
with higher Ro values, thereby indicating we should observe 
at least some evidence of a correlation. 

There appears to be no emergent correlation between Ro 
and radius deviations, at least for the present sample of data. 
We find no significant difference in the distribution of data 
points observed as a function of Ro. However, inspection of 
Figure [3] highlights the need for more data in order to draw a 
definitive conclusion. Selection of Ro = 0. 1 has the unfortu- 
nate effect of creating a bin with a small sample population, 
potentially affecting the statistics. Field stars in wide binary 
systems are a good starting point for studies interested in pop- 
ulating the high Ro region of Figure [3] Caution must further 
be taken as secondary stars in the longer period (> 15d) sys- 
tems in our sample do not have independently measured rota- 
tion periods, meaning they could potentially have different Ro 
values than are presented here. Visual inspection of Figure[3] 
leads us to the same conclusion, if we ignore the two most 
discrepant points. Since we were comparing MAE values, the 
outliers do not significantly affect the results of the statistical 
analysis. 



For fully convective stars, an a-co dynamo cannot operate due to the 
lack of a tachocline. Instead, it is thought that an a- dynamo can efficiently 
generate a magnetic field. Since it is the a mechanism that is related to the 
strength of the Coriolis force, the Rossby number should be just as applicable 
to fully convective stars iChabrier & Kiikei 2006; Browning 2008). 
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Figure 4. Relative error between observational and model radii for stars with detected x-ray emission. X-ray data is drawn from the ROSAT All-Sky Survey and 
combined with the radius residuals derived from this study. Data are shown as maroon filled circles. Illustrated are two least-square regressions performed on the 
data. The light-blue, dashed line demonstrates a non-negligible slope of ~ 25 ± 17 across all of the data, while the indigo, dash-dotted line excludes the two most 
discrepant points, UV Psc B and YY Gem. 

Arguably, the MAE is not an effective measure of the de- ^Wl^ 
gree of inflation of each sample as it treats deflated radii the 
same as inflated radii. Instead, the actual direct mean may 
provide a more compelling statistical measure. Thus, we ran 
our statistical analysis on the direct mean error. Overall, the 
typical degree of inflation among the radii of low-mass stars 
was found to be about 1.6%. No statistically significant cor- 
relations with either Porb or Ro were uncovered. The most 
significant result was found for the period cut at 1.5d, where 
the K-S test indicated a significant difference. However, the 
MC method produced a difference in the mean error of the 
two populations at 2.2o, below our significance threshold of 
3.0o. All other bin divisions yielded results significant at only 
about the la level. 

Curiously, if we consider only data points with measured 
x-ray fluxes (see below), there is evidence that systems with 
lower Ro values may have larger radius discrepancies. How- 
ever, we are then prompted to explain why the trend does not 
continue to higher Ro values, a question we are currently not 
equipped to answer There is still some ambiguity with the 
presence of the points near Ro = 0.01, which appear to con- 
tradict the presence of any definite correlation. Deeper x-ray 
observations of x-ray faint low-mass DEBs will clarify this 
ambiguity. 



iLopez-MoralesI (120071) previously performed such a com- 
parison and found a clear correlation between and radius 
deviations. Her comparison was performed under different 
modeling assumptions, which may have influenced the re- 
sults. Specifically, she compared all observations to a 1 Gyr, 
solar metallicity isochrone from Baraffe et al. (1998). With- 
out variation in age and metallicity, the observed radius dis- 
crepancies may have been over estimated (or under estimated 
in the case of YY Gem). 

A subsample of sixteen DEBs from this study have identi- 
fied x-ray counterparts in the ROSAT All-Sky Survey Bright 
and Faint Source Catalogues ( Voaes et al. 1999, 2000). Our 
analysis follows that of Lopez-Morales (2007) whose previ- 
ous analysis contained a fraction of our current x-ray detected 
sample. X-ray count rates were c onverted to x- r ay flu xes ac- 
cording to the formula derived bv lSchmitt et alj (Il995h . 



(5.30HR + 8.31) x lO^'^X, 



(5) 



6.2.2. X-ray Activity 

An interesting extension of our discussion in the previous 
subsection, is to compare our derived radius deviations with 
the observed x-ray to bolometric luminosity ratio (hereafter 
R^ = Lx/Lhoi). Since no correlation was noted with Ro, we 
expect that no correlation will be observed with R^, a s it has 
been shown to be intimately connected with Ro (.Wright et al.l 



where HR is the x-ray hardness ratiqj, Xa- is the x-ray count 
rate, and F^- is the x-ray flux. Luminosities in the x-ray spec- 
trum were computed with distances determine d from either 
Hippa rcos parallaxes (prefeiTed when available: lvan LeeuwenI 
l2007h or near-infrared photometry. 

Photometric distances were estimated using the luminos- 
ity calculated from the Stefan-Boltzmann relation assuming 
the observationally determined radius and effective temper- 
ature. Absolute magnitudes were then derived using the 
PHOENIX AMES-COND model atmospheres, adopting the 

There are typically two hardness ratios listed in the ROSAT catalog, HRl 
and HR2. The Schmitt et al. 1 1995) formula requires the use of HRl. 
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best fit isochrone metallicities. In an effort to reduce er- 
rors introduced by the theoretical atmospheres, distances were 
computed using the average distance modulus derived from 
2MASS J and K band photometry (Cutri et al. 2003). 

We derived values for all sixteen stars in our x-ray 
sample to ensure internal consistency. Slig ht discrepancies 
betwe en values presented here and those of iLopez-MoralesI 
(|2007') are due entirely to differences in the adopted distances. 
Attributing the x-ray flux contribution from each star in an 
DEB system is a difficult task. As such, Lopez-Morales per- 
formed her analysis using three reported emp irical scaling re- 
lations (iPallavicini etaLiri981l:lFleming et al.lll989h : 

• Case 1 - each component contributes equal weight. 

• Case 2 - proportional to the respective rotational veloc- 
ity, V sin /, of each component. 

• Case 3 - proportional to the square of the rotational ve- 
locity, (vsin/)^, of each component. 

We present results from all three cases in Figured Imme- 
diately we notice that the size of the stellar radius deviations 
appears to correlate with R^^. A linear least-square regression 
performed on each data set (light-blue dashed line in Figure|4]l 
suggests the slope for each case (1 - 3) is 26, 24, and 22 (± 
17), respectively, all with reduced-x^ values of ~8. Pearson 
r coefficients were 0.72, 0.69, and 0.63, respectively. The 
statistics suggest that the likelihood of uncorrelated sets of 
data producing these particular correlations are 0.2%, 0.3%, 
and 0.9%, for case 1, 2, and 3, respectively. We can therefore 
rule out the null hypothesis with greater than 99% confidence. 

Visually, however, we notice that the correlation is largely 
driven by the presence of YY Gem and UV Psc B, located in 
the upper-right region of each panel in Figure|4] If we were to 
remove those three points (YY Gem appears as a single point), 
the correlation vanishes and we only observe an offset from 
the zero point (Indigo dash-dotted line in Figure |4]i. Further- 
more, the distance to YY Gem is highly uncertain. Assuming 
it is associated with Castor AB provides a distance of about 
15 pc (see AppendixiBli. but our photometric analysis, as de- 
scribed above, places YY Gem at a distance of approximately 
1 1 pc, 4 pc closer to the Sun than the Castor AB system. In- 
stead of selecting a single distance estimate, we averaged the 
two estimates and adopted d — 13 ± 2 pc, wh i ch als o happens 
to be the distance assumed by Delfosse et al.l (120031) . 

The fact that our statistical correlation critically hinges 
upon three points, two of which are strongly distance de- 
pendent, is a cause for concern and implies that the statis- 
tics should be interpreted with care. If we believe the strong 
statistical correlation, then we are presented with a situation 
where the rotational data and the x-ray data disagree. This 
may be due to the physics underlying the stellar dynamo or 
those underlying x-ray saturation. However, there are two 
further interpretations that are contingent upon the role sys- 
tematics play. 

First, we have that the correlation is entirely real and the 
presence of non-inflated stars with Rx ^ 0.0007 are outliers 
in the relation. Second, systematics play a large role, as pro- 
posed in Section |5] Here, the truly deviant stars exhibit very 
strong x-ray emission {R^ > 0.001), while non-inflated stars 
show lower, varying levels of x-ray emission. For this view 
to hold, the relation between the level of radius inflation and 
magnetic field strength can not be linear. Strong magnetic 



fields would induce significant radius inflation while moder- 
ate and weak fields would produce little or no inflation. 

Accepting, on the other hand, that the statistical correlation 
is spurious, we are left with a picture that is coherent with our 
rotation analysis. Namely, magnetic activity may not be the 
leading cause for all of the observed inflated radii. Here, sys- 
tematics may still play a role in producing stars that appear 
inflated, but that are consistent with the models, leaving a few 
discrepant stars. Naturally inefficient convection, potentially 
dependent on particular stellar properties, may be operating. 
However, magnetic activity cannot be fully ruled out, as we 
do not yet have a fully self-consistent description of the in- 
teraction of magnetic fields with the stellar interior and atmo- 
sphere. It may be that magnetic fields acting within YY Gem 
and UV Psc have a more noticeable affect on stellar structure, 
as higher mass stars are affected more by changes to the prop- 
erties of convection (see Appendix lAli. 

We tend to favor a hybrid interpretation. Here, most of the 
observed "inflation" is an artifact of unaccounted-for system- 
atics, but significantly discrepant stars (YY Gem, UV Psc) 
are associated with very strong magnetic activity. We be- 
lieve that CM Dra probably fits into the latter category due 
to a push in the literature towards subsolar metallicity. Effects 
of a large-scale magnetic field are presumably mass depen- 
dent, with higher mass stars showing a greater propensity to 
become inflated owing to their sensitivity to changes in con- 
vective properties. Whether there is a characteristic magnetic 
field strength that induces substantial radius deviations is un- 
clear Dynamo saturation and the saturation of magnetic ac- 
tivity, as evidenced by the flattening of the Ro-Rj relation in 
Wrig ht et al.l (1201 ll) . is not yet fully understood, but may pro- 
vide further insight into the apparent disagreement between 
our rotation and x-ray analyses. 

Finally, clarity will be obtained with a better distance mea- 
surement to YY Gem, either from ground-based parallax pro- 
grams or with eventual results from Gaia. Of all the points in 
FigurelH the existence of a positive correlation is most depen- 
dent upon the distance to YY Gem. An accurate distance, as 
well as a reanalysis of its association with the Castor system, 
has the ability to not only relieve the ambiguity present in the 
x-ray data, but also to provide insight into the necessary con- 
straints for the system (i.e., is YY Gem really about 400 Myr 
old?). As further low-mass DEB systems are discovered, x- 
ray observations are strongly encouraged in order to develop a 
coherent picture of how radius deviations correlate with mag- 
netic activity. 

7. SUMMARY 

This study focused on reevaluating the current state of 
agreement between the theoretical and observational low- 
mass, main sequence MR relationship. The DEB systems 
used in the analysis were required to have quoted random un- 
certainties in the mass and radius below 3% in order to pro- 
vide an effective test of stellar evolution models. A large grid 
of DSEP models were computed with variation in age and 
metallicity characteristic of the local galactic neighborhood. 
Best fit isochrones were derived by allowing the age and 
metallicity to be optimized while maintaining the constraint 
that the system be coeval with a single composition. DEBs 
with reliably determined ages or metallicities were compared 
with a restricted set of isochrones in the range allowed by the 
observational priors. 

Overall, we find 92% of stars in our sample are less than 
4% discrepant with the models, largely representing a factor 
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of two improvement over the canonical 5 - 15% deviations. 
Our results suggest that low-mass stars with radii that deviate 
significantly from model predictions are exceptions to gen- 
eral agreement. Discrepancies may also be the result of unac- 
counted for systematic uncertainties (i.e., starspots) that may 
as large as 4%. With uncertainties as large as the typical ra- 
dius deviations, we find it difficult to draw the firm conclusion 
suggesting that models are in broad disagreement with obser- 
vations. Instead, we are left with a situation where the obser- 
vational uncertainties may be too large to provide an adequate 
test of stellar models. The combination of random and sys- 
tematic uncertainties for the sample of low-mass DEBs must 
be constrained and minimized below the 2% level before ac- 
curate model comparisons may be made. 

Radius correlations with orbital (rotational) period, Ro, and 
Rx were also considered. No distinct trends were identified 
with either orbital period or Ro. However, we find evidence 
for a strong corr elation between radius d eviations and Rx (pre- 
viously noted b v iLopez-MoralesI (l2007h ) in contradiction with 
ou r Ro analysis. The tre nd is not as tight as that derived 
by iLopez-MoralesI (120071) . owing to the age and metallicity 
variations allowed by our analysis. This correlation is also 
largely contingent upon the veracity of the distance estimate 
to YY Gem. Accurately determining the distance to YY Gem 
and evaluating its association with the Castor AB quadruple 
would alleviate much of the uncertainty. 

Finally, we must not leave the theoreticians out of the spot- 



light. The degree to which a magnetic field can alter the in- 
terior structure of low-mass stars is still only partially known 
and further investigations are required. Development of mod- 
els with self-consistent magnetic field perturbations will begin 
to shed light on this unknown. Comparisons between pre- 
dicted magnetic field strengths from self-consistent models 
and observational data (either direct or indirect) will provide 
a measure of the validity of the ability of magnetic fields to 
inflate stellar radii. Whatever the final solution may be to this 
long-standing problem, it is now apparent that the level of in- 
flation required by theory is not a ubiquitous 5 - 15%, but 
only so in extreme cases. 

We are grateful to G. Torres and J. Orosz for insightful dis- 
cussions and access to preliminary data results. The authors 
would also like to extend gratitude to A. Irwin for his work on 
the open-source project FreeEOS. GAF wishes to thank the 
William H. Neukom 1964 Institute for Computational Science 
for their generous support. BC and GAF would both also like 
to acknowledge the support of the National Science Founda- 
tion (NSF) grant AST-0908345. This reseai'ch has made use 
of NASA's Astrophysics Data System, the SIMBAD database, 
operated at CDS, Strasbourg, France, and the ROSAT data 
archive tools hosted by the High Energy Astrophysics Science 
Archive Research Center (HEASARC) at NASA's Goddard 
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APPENDIX 

A. VARIABLE MIXING-LENGTH MODELS 
A.l. Isochrone Grid 

A second set of isochrones was generated in order to address the idea that low-mass stars may possess inflated radii due 
to inefficient convective energy transport. It has been posited that convection within low-mass stars may either be inherently 
inefficient or that other physical processes (i.e., magnetic fields) act to reduce the ability of convection to effectively transport 
energy. We do not attempt to prescribe a physical mechanism associated with this conjecture. Instead, we attempt to parametrize 
convection in such a way so as to reduce the efficiency of convection in the lower mass regimes while still maintaining our solar 
calibration, necessary to properly model the Sun. 

The second grid of isochrones was computed with a mass-dependent mixing-length, henceforth referred to as "variable OLmlt" 
models. Convection for these models was parametrized with a smooth quadratic function of the form 

a ]^a[ +b. (Al) 

Selection of a quadratic was arbitrary and carries no physical justification, except to produce low-mass stars with relatively 
inefficient convection compared to those in the standard model case. 

Coefficients were determined by matching the convective mixing-length to predetermined values at two different masses. 
Since the overall structure of very low-mass stars is rather insensitive to the precise value of the mixing-length, we anchored 
<^MLT = 1 00 at M = 0. IM©. The other end of the mass spectrum is constrained by our need to satisfy our solar calibration. Thus, 
at M = 1 .OMq, the convective mixing-length was fixed to Umlt = 1 -94. Subsequently, 

af—^ =0.949 f—^ +0.991. (A2) 

A relative comparison between models computed with a solar calibrated mixing-length and those generated with a parametrized 
mixing-length is presented in Figure |5]A. 

Our isochrone grid for these models covered a fraction of the parameter space compared to the standard model set. There was 
a total of 12 isochrones for the variable Umlt models: [Fe/H] = {-0.5, 0.0, H-0.2} dex with ages = {1.0, 3.0, 5.0, 8.0} Gyr. Below 
the fully convective boundary, the models are rather insensitive to the mixing-length, as stated above. For higher masses, the 
reduced mixing-length increases the model radii by up to 3%. As designed, a Sun-like star is unaffected by this mixing-length 
prescription. We also see that above IMp, the model radii begin to decrease due to an increased mixing-length as a result of our 
parametrization. 
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Table 5 

Best fit isochrone with a mass-dependent convective mixing-length. 



Quoted Fixed 3% 



Star Name 


Age 


[Fe/H] 


bR/Robs 




Fit 


Age 


[Fe/H] 






Fit 


UV Psc A 


8.0 


0.00 


-0.0072 


-0.349 


No 


8.0 


0.00 


-0.0072 


-0.241 


No 


UV Psc B 






0.0866 


4.019 








0.0866 


2.888 




IM Vir A 


8.0 


+0.20 


-0.0032 


-0.209 


No 


8.0 


+0.20 


-0.0032 


-0.105 


Yes 


IM Vir B 






0.0238 


1.245 








0.0238 


0.792 




KID 6131659 A 


3.0 


-0.50 


-0.0113 


-5.876 


No 


3.0 


-0.50 


0.01 13 


-0.378 


Yes 


KID 6131659 B 






-0.0190 


-9.352 








-0.0190 


-0.634 




RX J0239.1-1028 A 


8.0 


-0.50 


-0.0075 


-1.394 


No 


8.0 


-0.50 


-0.0075 


-0.251 


Yes 


RX 10239.1-1028 B 






0.0039 


1.356 








0.0039 


0.129 




Kepler- 16 A 


1.0 


0.00 


-0.0085 


-4.227 


No 


1.0 


0.00 


-0.0085 


-0.282 


Yes 


Kepler- 16 B 






0.0244 


11.061 








0.0245 


0.815 




GU Boo A 


8.0 


-0.50 


0.0085 


0.333 


Yes 


8.0 


-0.50 


0.0085 


0.283 


Yes 


GU Boo B 






0.0190 


0.740 








0.0190 


0.632 




YY Gem A 


0.4 


0.00 


0.0757 


8.226 


No 


0.4 


0.00 


0.0757 


2.523 


No 


YY Gem B 






0.0757 


8.226 








0.0757 


2.523 




MG 1-5 06664 A 


1.0 


0.00 


-0.0017 


-0.964 


No 


1.0 


0.00 


-0.0017 


-0.057 


Yes 


MG 1-5 06664 B 






-0.0099 


-5.065 








-0.0099 


-0.329 




MGl-1 16309 A 


8.0 


+0.20 


-0.0126 


-1.743 


No 


8.0 


+0.20 


-0.0126 


-0.421 


Yes 


MGl-1 16309 B 






0.01 15 


1.534 








0.0115 


0.384 




MGl-1819499 A 


5.0 


-0.50 


0.0240 


6.839 


No 


5.0 


0.00 


0.0413 


1.376 


No 


MGl-1819499 B 






-0.0575 


-9.587 








-0.0396 


-1.320 




NSVS 01031772 A 


8.0 


+0.20 


-0.0121 


-2.276 


No 


8.0 


+0.20 


-0.0121 


-0.404 


No 


NSVS 01031772 B 






0.0352 


5.965 








0.0352 


1.172 




MGl-78457 A 


3.0 


0.00 


-0.0023 


-0.143 


Yes 


3.0 


0.00 


-0.0023 


-0.075 


Yes 


MGl-78457 B 






0.0005 


0.025 








0.0005 


0.016 




MG 1-646680 A 


1.0 


0.00 


-0.0319 


-2.428 


No 


1.0 


0.00 


-0.0319 


-1.063 


No 


MG 1-646680 B 






0.0159 


1.129 








0.0159 


0.529 




MGl-2056316 A 


1.0 


+0.20 


-0.0083 


-1.829 


No 


1.0 


+0.20 


-0.0083 


-0.277 


Yes 


MG1-2056316B 






0.0112 


2.095 








0.0112 


0.373 




CU Cnc A 


8.0 


+0.20 


0.0127 


0.998 


Yes 


8.0 


+0.20 


0.0127 


0.423 


Yes 


CU Cnc B 






-0.0028 


-0.116 








-0.0028 


-0.092 




LSPM 111 12+7626 A 


8.0 


+0.20 


-0.0038 


-0.278 


No 


8.0 


+0.20 


-0.0038 


-0.125 


Yes 


LSPM 111 12+7626 B 






0.0255 


1.652 








0.0255 


0.850 




KOI- 126 B 


3.0 


+0.20 


0.0020 


0.365 


Yes 


5.0 


0.00 


0.0010 


0.033 


Yes 


KOI-126C 






-0.0017 


-0.307 








-0.0014 


-0.047 




CM Dra A 


5.0 


0.00 


0.0282 


3.757 


No 


5.0 


0.00 


0.0282 


0.939 


No 


CM Dra B 






0.0333 


4.441 








0.0333 


1.111 





A.2. Results 

Non-standard, variable mixing-length models lead to only slightly better results over the standard models when considering 
the quoted random uncertainties (see Table|5]). We again find that radius deviations are largely reduced to below 4%, seen clearly 
in Figure |5]B. The residuals do appear to be more tightly clustered around the zero point, although only slightly, suggesting an 
overall better agreement between models and observations. This is reinforced by a slightly lower MAE of 2. 1 % across the entire 
sample. 

There was evidence in Figure|5]B of a possible tendency for the best fit theoretical isochrone to over predict the observed radius 
at higher masses and under predict the radius for lower mass stars. A least squares regression on the data finds a linear slope of 
—0.031 ± 0.013, significant at the 2.5a level (not shown). Specifically, half of the DEB systems that were not formally fit are 
characterized by an isochrone that over-predicts the higher mass primary but under-predicts the lower mass secondary. This is 
only an artifact of our functional parametrization of Umlt, suggesting convection was too heavily suppressed at higher masses 
and under-suppressed in the low-mass regime. 

We subjected the resulting data to the same statistical analyses described in Section 15.11 We again compared the radius 
deviations to the rotational period and Rossby number (Figure |5]C and D). The strongest hint of a correlation was found when 
applying a cut at a period of 1.5d, although it was only significant at the 2.7a level, below our 3a significance threshold. 
Performing the period cut at l.Od and 2.0d yielded results that were even less significant (2.4a and 1.9a, respectively), providing 
the same qualitative result as standard model scenario. Similarly, there is no strong correlation observed between the radius 
deviations and the Rossby number. 

Again, artificially fixing the uncertainties at 3% tends to add ambiguity as to whether there are any real discrepancies for 
a majority of the systems. As was suggested throughout our study of the standard models, systematic uncertainties must be 
minimized to below 2% before an accurate comparison between models and observations may occur. 

B. AGE AND METALLICITY PRIORS 

Mentioned in Section 14.41 it was determined that four of the DEBs in our sample had age and metallicity priors reliably 
determined. These four were a subset of a total of eight DEB systems that had either an estimated age or metallicity in the 
literature. Below, we provide our analysis and reasoning for our acceptance or rejection of each quoted prior 
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Figure 5. Top-left: Radius variations induced by our prescription of a mass-dependent cImlt (Eq.|A2j. Sliown are 1 Gyr (blue - solid) and 8 Gyr (maroon - 
dashed) isochrones computed with a solar heavy element composition (GS98). Positive values indicate that the modified oImlt models are inflated compared to 
the sol ar calibrated models. The "blitf' observed near the fully convective boundary is most likely due to a ^^He instability described by van Saders & PinsonneauIJ 
120121) . Top-right: Same as Figure|T|A excepted for the modified oLmlt models. Bottom-left: Same as Figure|2]B, although we only present the short period 
systems. The distribution of long period systems is similar to Figure|2]A. Bottom-right: Identical to Figure|3] 

B.l. UVPsc 

The age of the UV Psc system is quoted by the TAGIO review to be 7.9 Gyr Popper's original paper describing the character- 
istics of the UV Psc system does not provide any evidence to support an age estimate (Popper 1997). An age of approximately 8 
Gyr is typically assigned to the system due to the fact that stellar evolution models predict the physical properties of the primary 
star, UV Psc A, at that age. However, there has yet to be any set of models that can place both components on a consistent, coeval 
isochrone. Since the estimated age of UV Psc A appears to be derived from stellar evolution models, we allowed our models to 
independently determine the most acceptable age. 

B.2. IMVir 

[Morales et all ( l2009bl) attempt to determine the metallicty of IM Vir by applying various photometric metallicity relations. The 
result of their efforts was that they found all of the various empirical methods quote different values with large uncertainties. 
Values for [Fe/H] vary from the metal -poor end with [Fe/H] = -0.8 up to a super-solar value of [Fe/H] = +0.15. This range also 
fits nicely within the set of model metallicity values selected for this study. Since the cited metallicity range would not provide 
any additional constraint on our analysis, we rejected the metallicity prior. 

B.3. Kepler-16 

Kepler- 16, the first confirm ed binary system with a circumbinary planet, was provided with a metallicity estimate in its discov- 
ery paper (iDovle et al.ll20Tlh . A metallicity of [Fe/H] = -0.3±0.2 was determined spectroscopically. The authors indicated the 
spectroscopic analysis was performed on the K-dwarf primary and that the general method was similar to that applied to KOI- 126 
using Spectroscopy Made Easy (SME; Valenti & Piskunov 1996; Carter et al. 201 1). Reliability is lent to the method, in general, 
due to its success at deriving the metallicity of KOI- 126 A. We therefore accepted the quoted metallicity prior 

B.4. GUBoo 

The age estimate provided by 'Lopez-M orales & RibasI (12005 b in their characterization of GU Boo was primarily based on 
kinematics. Specifically, they conclude that GU Boo is an isolated system and the vertical component of its motion provides 
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a hint that it has undergone perturbations due to disk heating processes. Assuming that the system has been subjected to disk 
heating, one can only infer that the system has an age gr eater then 10^ yr , the typical timescale for dynamical perturbations 
associated with an objects orbit around the galactic center (!Soderblomll2010l) . Unfortunately, no further constraints were able to 
be placed on the age of the system. However, it is not possible to rule out the scenario that GU Boo was dynamically ejected from 
its stellar nursery. It would then appear that no reliable age estimation exists, prompting us to reject the age prior for this system. 

B.5. YYGem 

YY Gem is thought to be physically associated with the Castor AB (a Gem) quadruple system. All three systems were 
found to be gravi tationally bound base d on a statistical analysis performed using a three-body interaction code dAnosova et alj 
119891) . While the lAnosova et all (|1989[) analysis was performed with pm-Hipparcos proper motions as their initial conditions, it 
is unlikely that the results will be effected at the level necessary to unbind the systems. 

Castor A and B are themselves both binaries. The primary in both systems is an A star and both are thought to have an M dwarf 
companion. Therefore, physical properties derived from spectroscopic and photometric observations of the two A stars will be 
essentially unaffected by the presence of their companions. Placing the two primaries on an Mv-Log(refl) H-R diag ram , we u sed 
DSEP to derive an age of about 400 Myr This age is consistent with the average age of 370 Myr derived by iTorres & RibasI 
(l2002i) . who modeled the A-star primaries using multiple stellar evolution codes. 

The A stars also lend t hemse l ves w ell to a spectroscopic determination of the metallicity. Unfortunately, there appears to be 
only one result reported. Smith' (1914) estimates an average metallicity for the two A stars to be about [Fe/H] = +0.1 relative to 
Vega. As described in Torres & Ribas ( 2002), this implies a rather uncertain metallicity relative to the Sun of [Fe/H] = H-0. 1 ± 0.2. 
Despite this, we adopt this metallicity constraint due to its rather large uncertainty, which should presumably encompass the true 
value. 

B.6. CUCnc 

When CU Cnc was originally investigated by I RibasI ( l2003h . it was found to have a space motion very similar to that of the 
Castor sextuple system (a Gem, YY Gem) and, subsequently, the proposed "Castor moving group." For this reason, CU Cnc was 
deemed to be associated with the Castor moving group implying an age and metallicity similar to the Castor sextuple. However, 
determining membership of a moving group is complicated and has often lead to ambiguous results concerning the coeval nature 
of the group. It is not uncommon for members of the same kinematic moving group to have different metallicities, implying that 
members of a defined moving group may not have been born in the same galactic environment and, as such, are not necessarily 
coeval. The lack of an age estimation beyond its kinematic similarity to the Castor system led us to reject the age prior of CU 
Cnc. 

B.7. KOI-126 

KOl-126 is a hierarchical triple eclipsing binary recently whose discovery was recently announced (ICarter et al .ll201 ll). The 
quoted metallicity prior of [Fe/H] = H-0. 15 ± 0.08 was determined using SME, as was mentioned above in the discussion of Kepler 
16. Assuming this metallicity allowed for a relatively precise age constraint (~'4±1 Gyr) to be placed on the primary, a 1.35 Mq 
subgiant. Combini ng both the age and metallicity information led to two low-mass companions to also be fit by stellar models 
dFeiden et alj|201 ll) . Thus, there appears to be little question about the validity of the age and metallicity estimations, leading us 
to adopt the given priors for our study. 

B.8. CMDra 

Finally, we consider the very well studied CM Dra system. Spectroscopic observations of the system have produced varying 
results for the metallicity of th e system, but all appear to be consistent with -1 < [Fe/H] < ('Viti et al. 1997, 2002; Morales et alj 
l2009atlKuznetsov et alj|2012h . Due to the difficulties in modeling, and thus fitting, the entire SED of an M dwarf, we determined 
there was no particular reason to strongly favor one metallicity result over another. 

An age was determined for the system through th e use of white d warf (WD) cooling tracks combined with stellar evolu- 
tion models o f the a pproximate WD progenitor (Mor ales et al.ll20()9ah . The WD cooling age was found to be 2.38±0.37 Gyr 
[Morales et alJ (12009a) predicted a mass for the progenitor star of M = 2.1 ± QAMp, and derived a total (stellar model + WD 
cooling) age of 4.1 ±0.8 Gyr. For consistency, we calculated the lifetime of the progenitor star using DSEP with the full available 
suite of physics. Accounting for the metallicity constraint defined above, we derived an age of 3.7±1.2 Gyr and, in doing so, 
found no reason to reject, or modify, the age prior 

REFERENCES 



Andersen, J. 1991, A&A Rev., 3, 91 

Anosova, J. R, Orlov, V. V., & Chemyshev, M. V. 1989, Soviet Ast. Let., 15, 
237 

Bahcall, J. N., Basu, S., Pinsonneault, M., & Serenelli, A. M. 2005, ApJ, 
618, 1049 

Baraffe, 1., Chabrier, G., Allard, R, & Hauschildt, R H. 1997, A&A, 327, 
1054 

Baraffe, 1., Chabrier, G., Allard, E, & Hauschildt, R H. 1998, A&A, 412, 
403 

Barnes, S. A. 2010, ApJ, 722, 222 



Berger, D. H., Gies, D. R., McAlister, H. A., et al. 2006, ApJ, 644, 475 
Bjork, S. R. & Chaboyer, B. 2006, ApJ, 641, 1102 
Browning, M. 2008, ApJ, 676, 1262 

Burrows, A., Heng, K., & Nampaisarn, T. 201 1, ApJ, 736, 47 

Carter, J. A., Fabrycky, D. C, Ragozzine, D., et al. 2011, Science, 331, 562 

Chaboyer, B., Fenton, W. H., Nelan, J. E., Patnaude, D. J., & Simon E E. 

2001, ApJ, 562, 521 
Chaboyer, B. & Kim, Y.-C. 1995, ApJ, 454, 767 
Chabrier, G. & Baraffe, 1. 1997, A&A, 327, 1039 
Chabrier, G. & Baraffe, 1. 2000, ARA&A, 38, 337 



16 



Feiden & Chaboyer 



Chabrier, G. & Kiiker, M. 2006, A&A, 446, 1027 
Chabrier, G., Gallardo, J., & Baraffe, I. 2007 A&A,472, L17 
Chandrasekhar, S. 1933. MNRAS, 93, 390 

Cutri, R. M.. Skrutskie, M. R, van Dyk, S., et al. 2003, VizieR Online Data 
Catalog, 2246, 

Delfosse, X., Forveille, T., Segransan, D., et al. 2000, A&A, 364, 217 
Dotter, A.. Chaboyer, B., Jevremovic, D., et al. 2007, AJ, 134, 376 
Dotter, A., Chaboyer, B., Jevremovic, D., et al. 2008, ApJS, 178, 89 
Doyle, L. R., Carter, J. A., Fabrycky, D. C, et al. 2011, Science, 333, 1602 
Feiden, G. A., Chaboyer, B., & Dotter, A. 2011, ApJ, 740, L25 
Ferguson, J. W., Alexander, D. R., Allard, F., et al. 2005. ApJ, 623, 585 
Fleming, T. A., Gioia, I. M., & Maccacaro, T. 1989, ApJ, 340, 101 1 
Gough, D. O. & Tayler, R. J. 1966, MNRAS, 133, 85 
Grevesse. N., & Suaval, A. J. 1998, Space Sci. Rev, 85, 161 
Guenther, D. B., Demarque, R, Kim, Y.-C, & Pinsonneault, M. H. 1992, 
ApJ, 387, 372 

Hauschildt, R H., Allard, F, & Baron, E. 1999a, ApJ, 512, 377 
Hauschildt, P. H., Allard, F, Ferguson, J., Baron, E., & Alexander, D. R. 

1999b, ApJ, 525, 871 
Hoxie, D. T. 1970, ApJ, 161, 1083 
Hoxie, D. T. 1973, A&A, 26, 437 
Iglesias, C. A. & Rogers, F J. 1996, ApJ, 464, 943 
Irwin, J. M., Quinn, S. N., Berta, Z. K., et al. 201 1, ApJ, 742, 123 
Jackson, R. J., Jeffries, R. D., & Maxted R R L. 2009, MNRAS, 399, L89 
Kim, Y.-C. & Demarque, R 1996, ApJ, 457, 340 

Kjurkchieva, D. R, Marchev, D. V., Heckert, R A., & Ordway, J. I. 2005, AJ, 
129, 1084 

Korn, A. J., Grundahl, R, Richard, O., et al. 2007, ApJ, 671, 402 

Kraus, A. L., Tucker, R. A., Thompson, M. 1., Craine, E. R., & Hillenbrand, 

L. A. 2011, ApJ, 728, 48 
Kuznetsov, M. K., Pavlenko, Y. V., Jones, H. R. A., & Pinfield, D. D. 2012, 

Advances in Astronomy and Space Physics, 2, 15 
Lopez-Morales, M. 2007, ApJ, 660, 732 
Lopez-Morales, M., Orosz, J. A., Shaw, J. S., et al. 2006, 

arXiv:astro-ph/0610225 
Lopez-Morales, M. & Ribas, I. 2005, ApJ, 631, 1120 

Lopez-Morales, M. & Shaw, J. S. 2007, in ASP Conf. Ser. 362, The Seventh 
Pacific Rim Conference on Stellar Astrophysics, ed. Y. W. Kang et al. 
(San Prancisco, CA: ASP), 26 

MacDonald, J. & MuUan, D. J. 2012, IVINRAS, 421, 3084 



Morales, J. C, Gallardo, J., Ribas, 1., et al. 2010, ApJ, 718, 502 
Morales, J. C, Ribas, I., & Jordi, C. 2008, A&A, 478, 507 
Morales, J. C, Ribas, I., Jordi, C, et al. 2009, ApJ, 691, 1400 
Morales, J. C, Torres, G., Marschall, L. A., & Brehm, W. 2009, ApJ, 707, 
671 

Mullan, D. J. & MacDonald, J. 2001, ApJ, 559, 353 

Noyes, R. W., Hartmann, L. W., Baliunas, S. L., Duncan, D. K., & Vaughan, 

A. H. 1984, ApJ, 279, 763 
Pallavicini, R., Golub, L., Rosner, R., et al. 1981, ApJ, 248, 279 
Parker, E. N. 1975, ApJ, 198, 205 

Parker, E. N. 1979, Cosmical Magnetic Rields (Oxford: Clarendon Press) 

Popper, D. M. 1997, AJ, 114, 1195 

Reiners, A. 2012, Living Reviews in Solar Physics, 9, 1 

Reiners, A., Joshi, N., & Goldman, B. 2012, AJ, 143, 93 

Ribas, I. 2003, A&A, 398, 239 

Ribas, I. 2006, Ap&SS, 304, 89 

Richard, O., Michaud, G., & Richer, J. 2005, ApJ, 619, 538 
Schmitt, J. H. M. M., Pleming, T. A., & Giampapa, M. S. 1995, ApJ, 450, 
392 

Skumanich, A. 1972, ApJ, 171, 565 

Smith, M. A. 1974, ApJ, 189, 101 

Soderblom, D. R. ARA&A, 48, 581 

Spada, R & Demarque, P 2012, MNRAS, 422, 2255 

Thoul, A. A., Bahcall, J. N., & Loeb, A. 1994, ApJ, 421, 828 

Torres, G., Andersen, J., & Gimenez, A. 2010, A&A Rev., 18, 67 

Torres, G. & Ribas, 1. 2002, ApJ, 567, 1 140 

Vaiana, G. S., Cassinelli, J. P, Rabbiano, G., et al. 1981, ApJ, 244, 163 
Valenti, J. A. & Piskunov, N. 1996, A&AS, 118, 595 
van Leeuwen, R 2007, A&A, 474, 653 

van Saders, J. L., & Pinsonneault, M. H. 2012, |arXiv:1201.5381 l 
Viti, S., Jones, H. R. A., Maxted, P, & Tennyson, J. 2002, MNRAS, 329, 
290 

Viti, S., Jones, H. R. A., Schweitzer, A., et al. 1997, MNRAS, 291, 780 
Voges, W., Aschenbach, B., Boiler, Th., et al. 1999, A&A, 349, 389 
Voges, W., Aschenbach, B., Boiler, Th., et al. 2000, VizieR Online Data 
Catalog, 9029, 

von Braun, K., Boyajian, T. S., Kane, S. R., et al. 2012, arXiv: 1 202.00831 
Windmiller, G., Orosz, J. A., & Etzel, P B. 2010, ApJ, 712, 1003 
Wright, N. J., Drake, J. J., Mamajek, E. E., & Henry, G. W. 201 1, ApJ, 743, 
48 



